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Q ! ABSTRACT 

^ ' Context. The formation of massive stars is not well-understood and requires detailed observa- 

tional studies in order to discriminate between the different proposed star formation models. 
Aims. We have searched for massive molecular outflows in a sample of high-mass star forming re- 
gions, and we have characterised both the outflow properties and those of their associated molec- 



> 

00 . 

, ular clumps. With a sample composed largely of more luminous objects than previous ones, this 

work complements analogous surveys performed by other authors by adding the missing highest 
luminosity sources. 



' Methods. The sample under study has been selected so as to favour the earliest evolutionary 

' phases of star formation, and is composed of very luminous objects {L\,„\ > 2 x \(f Lq and up 

to ~10* Lq), possibly containing 0-type stars. Each source has been mapped in '^^CO(2-l) (an 



outflow tracer) and C'*'0(2-l) (an ambient gas tracer) with the lRAM-30m telescope on Pico 
Veleta (Spain). 

\ Results. The whole sample shows high- velocity wings in the ''C0(2-l) spectra, indicative of 

outflowing motions. In addition, we have obtained outflow maps in 9 of our 1 1 sources, which 
display well-defined blue and/or red lobes. For these sources, the outflow parameters have been 
derived from the line wing '^C0(2-l) emission. An estimate of the clump masses from the 
C'*0(2-l) emission is also provided and found to be comparable to the virial masses. From 
a comparison between our results and those found by other authors at lower masses, it is clear 
that the outflow mechanical force increases with the bolometric luminosity of the clump and with 
the ionising photon rate of the associated Hii regions, indicating that high-mass stars drive more 
powerful outflows. A tight correlation between outflow mass and clump mass is also found. 
Conclusions. Molecular outflows are found to be as common in massive star forming regions 
as in low-mass star forming regions. This, added to the detection of a few tentative large-scale 
rotating structures suggests that high-mass stars may generally form via accretion, as low-mass 
stars. 
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1. Introduction 

Massive stars have a high impact on their environments and play a major role in the evolution of 
the Galaxy. However, their formation mechanism is far from being well understood. There is a 
main theoretical problem that arises when studying massive star formation: unlike low-mass stars, 
stars with masses greater than ~8 Mq reach the zero-age main sequence (ZAMS) while they are 
still accreting material from the surrounding parental molecular cloud (Palla & Stabler [T993]l. The 
radiation pressure they develop at this point is expected to exceed the gravitational pressure of the 
collapsing material, thus halting the accretion process and preventing the star from gaining more 
mass. This leads to the conclusion that stars more massive than about 8 Mq cannot form, which 
is in clear disagreement with observations. The various models which have been proposed to solve 
this problem can be grouped into two main scenarios: accretion through massive circumstellar 
disks and/or sufficiently high accretion rates (e.g. Yorke & Sonnhalter l2002l Tan & McKee 120021 
Keto et al. I2003I I, and coalescence of lower-mass stars in very dense clusters (Bonnell, Bate & 
ZinneckerEOOll Bonnell & Bate|2002ll. 

Observing the early phases of the evolution of massive protostars is rendered more difficult 
by the high extinction, the relatively large distance (several kiloparsec) and by the considerable 
confusion. Massive stars do not form alone but in a cluster and the phenomena associated with 
the formation, say, of an O-star, have to be disentangled from those associated with the forma- 
tion of lower mass stars. Reviews of different aspects of the problem are given by Beuther et al. 
(I2007al l. Hoare et al. (l2007l l. McKee & Ostriker (|2007] |. and Zinnecker & Yorke (l2007l l. One gen- 
eral conclusion is that feedback due to outflows and ionising radiation from young massive stars 
may drastically influence the formation of still younger stars. Outflows and water or methanol 
masers are often moreover the most obvious signs of a young massive protostar It is therefore of 
considerable interest to determine the physical properties of the outflows detected in high-mass star 
forming regions (SFRs). 

Outflows in high-mass SFRs have had considerable study (see e.g. Arce et al. I2007I I. Two in- 
teresting examples are the flows associated with G24.78-h0.08 and G3 1.41 +0.31. These are two 
high-luminosity sources showing signposts of high-mass star formation, such as water masers and 
ultracompact (UC) Hn regions, which have been observed in a variety of molecular tracers, both 
with single-dish telescopes and interferometers (e.g. Cesaroni et al. II 9941 Olmi et al. II 9961 Furuya 
et al. 120021 Beltran et al. I2004l l. The observations reveal disk-outflow systems, and suggest that 
high-mass stars in these two regions form via accretion inside a pc-scale molecular clump. It is 
clear from the results of these particular studies that only the combination of low- and high-angular 
resolution observations can provide us with the necessary information on the structure and kine- 
matics of high-mass SFRs. Despite the great importance of these findings, the formation scenario 
depicted for G24 and G3 1 does not necessarily have to represent the general case. Therefore, more 
examples like these need to be found and investigated in detail in order to gradually build up a sta- 
tistically significant collection of well-known massive SFRs, and thus understand better the process 
of high-mass star formation. 



* Based on observations carried out with the IRAM 30-m telescope on Pico Veleta (Granada, Spain). IRAM 
is supported by INSU/CNRS (France), MPG (Germany), and IGN (Spain). 
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Outflow studies to date have mostly been carried out using '^CO or SiO. Using SiO is problem- 
atic because it seems likely that the SiO abundance varies widely between outflows and perhaps as 
a function of velocity within a single outflow. Using '^CO is also difficult because of the effects 
of optical depth as illustrated by the work of Choi et al. ( I1993I I. who found optical depths of the 
order of 10 in the CO(3-2) line towards a couple of flows in high mass star forming regions. They 
also found that the optical depth dropped off sharply in velocity channels well separated from the 
ambient value. High optical depth in '^CO has also to be contended with in the lower J transitions 
and has the consequence that outflow mass estimates based on the assumption of low optical depth 
in '^CO are underestimates. Partly for this reason, we have decided to estimate outflow mass and 
momentum in this article using the '^CO(2-l) line and correcting for possible finite optical depth 
in this transition using the wings of C"*0(2-l). While we expect this technique will yield a more 
accurate outflow mass estimate than that from '^CO, one misses in this way the high velocity CO 
wings (full width at zero power above 100 km s ' in some of the Choi et al. sample) and thus we 
underestimate the mechanical energy output. 

As a starting point for this study, we have conducted a small survey towards 1 1 high-mass 
star forming regions. This allows us to derive physical properties of a reasonable number of ob- 
jects and, perhaps most importantly, to select the best candidates for follow-up observations at 
higher spatial resolutions. A search for massive molecular outflows at millimetric wavelengths is 
well-suited for a single-dish project, since high-mass outflows are extended structures with sizes 
of ~1 pc (e.g. Shepherd & Church well 1 1 9961 Beuther et al. I2002bl l. i.e. a few half -power beam 
widths (HPBWs) at the typical distances of these sources (1-10 kpc). Moreover, the detection of a 
bipolar molecular outflow in a high-mass SFR already suggests the presence of an accretion disk 
and encourages further investigation at higher-angular resolution to search for direct signatures of 
rotating disks/toroids and more collimated jets closer to the central object. 

With this idea in mind, we have carried out a '^CO(2-l) and C'^0(2-l) survey with the 
IRAM 30-m telescope (Granada, Spain) towards a sample of selected high-mass young stellar 
objects (YSOs) in their earliest evolutionary stages, to search for massive molecular outflows. We 
describe the sample and the observations in Sect. 2. The results of our single-dish data are presented 
in Sect. 3. In Sect. 4 we derive the physical parameters of the outflows and the molecular clumps 
hosting them, in Sect. 5 we discuss the general results and comment on some particular sources, 
and finally. Sect. 6 summarises our results and conclusions. 

2. Observations and data reduction 

2. 1 . The sample 

The sample observed in the present study is composed of 1 1 high-mass SFRs chosen from the 
H2O maser surveys by Hofner & Churchwell (1996) and Forster & Caswell (1999), and from 
the CH3OH maser survey by Walsh et al. (1997). All the sources in our sample have bolometric 
luminosities greater than 2 x lO'* Lq (i.e. possibly with 0-type stars) and present compact or no 
free-free emission, i.e. they are associated with either UC Hii or no detected Hii regions. 

In addition, they are associated with at least one of the following: (i) H2O masers, (ii) CH3OH 
masers, (iii) compact (sub)mm continuum emission. All these criteria ensure the selection of high- 
mass YSOs in their earliest evolutionary phases. 
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Table 1. The sample: coordinates used for the observations, distance, bolometric luminosity, ion- 
ising photon rate, methanol maser detection and outflow detection 



Source a (2000) 5(2000) d d^^{' Ltoi* A^Ly CH3OH masers?* HjO masers? Outflow? 



(hms) ("'") (kpc) (kpc) 10^ L© (10^' s"') (Y/N) (Y/N) (Y/N) 



GlO.47+0.03 


18 08 38.22 


-19 51 49.7 


10.8" 


2.9 


140' 


28'' 


Y 


Y" 


Y 


GlO.62-0.38 


18 10 28.70 


-19 55 49.7 


6.5" 


2.4 


110*^ 


91' 


Y 


Y" 


Y 


G16.59-0.06 


18 21 09.16 


-14 31 48.8 


4.7'^ 


4.2 


2.3"^ 


0.26" 


Y 


Y" 


Y 


G19.61-0.23 


18 27 38.15 


-11 56 38.5 


12" 


4.9 


170^ 


49' 


N 


yb 


Y 


G23.44-0.18 


18 34 39.25 


-08 31 38.8 


7.8'' 


3.4 


57s 


<0.02" 


N 


yd 


7 


G28.87+0.06 


18 43 46.24 


-03 35 30.4 


8.5'' 


4.2 


29" 


0.13'" 


Y? 


yd 


Y 


G29.96-0.02 


18 46 03.96 


-02 39 21.5 


9= 


4.5 


190^ 


31'' 


Y 


yb 


7 


G35.20-0.74 


18 58 12.98 


+01 40 37.6 


2.3'' 


6.8 


2.^ 


0.02" 


Y? 


yd 


Y 


G43. 89-0.78 


19 14 26.16 


+09 22 34.0 


4.2" 


6.2 


4.0' 


9.1' 


Y? 


yb 


Y 


G48.61+0.02 


19 20 31.19 


+ 13 55 24.9 


11.8" 


8.9 


130' 


4.4° 




yd 


Y 


G75.78+0.34 


20 21 44.10 


+37 26 39.5 


4.1" 


8.5 


24i 


3.9' 


Y? 


yb 


Y 



t Galactocentric distance 

* Luminosities and ionising photon rates corrected for revised distances 

* Pestalozzi et al. ( 12005 b . Sources with "?" have only single-dish observations; "-" denotes no observa- 
tions towards that source. 

References: " Pandian et al. ( |2008] |. Hofner & Churchwell {T996l l, Codella et al. ( [T997l l. Forster & 
Caswell ( fT999t , "= Hanson et al. i2002t . ' Hatchell et al. ( [2000] >. s Walsh et al. (1997), " Fuller et al. ( [200ll . 
' Kurtz et al. ( [1994) , J Klaassen & Wilson l [2007] >. ''Cesaroni et al. ( fT994] >. ' Wood & Churchwell ( fT989) , 
Cesaroni (priv. comm.), " Codella (priv. comm.) 

In Table[T]we present the list of the observed sources, their coordinates, distances and bolomet- 
ric luminosities. The last three columns show, respectively, the rate of ionising photons from the 
associated ultracompact (UC) Hii regions, the methanol maser detection as reported by Pestalozzi 
et al. (2005), and the molecular outflow detection according to our observations (see Sect. 3.1). 
Note that all the sources are associated with UC Hii regions but G23.44, where no ionised emission 
has been detected so far. For this source, an upper limit for Nly is given, determined from 3cr. Both 
Lboi and A^Ly have been corrected for revised distances, based on the references in Table [T] 

2.2. IRAM 30-m observations 

The IRAM 30-m telescope at Pico Veleta (near Granada, Spain) was used to map the 1 1 sources 
of our sample in the '■'CO(2-l) and C'^0(2-l) lines. We have chosen the less abundant '^CO 
isotopologue instead of the more commonly used outflow tracer '^CO in order to decrease the 
effects due to confusion with other unrelated clouds along the line of sight, which is likely, given 
the galactic distribution of the objects (see also the discussion in Sect. 5.1). The observations were 
carried out on September 23 and 24, 2006. The average T^y^ during the observations was ~1000 K 
for the '^CO(2-l) and ~1400 K for the C''^0(2-l) line. The pointing was checked about every 
hour by observing nearby planets and it was found to be accurate to within 4". 
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The maps have sizes of 4'x4' and were obtained with the 3x3 pixel receiver array HERA in the 
on-the-fly mode, and the autocorrelator VESPA. At the frequency of the lines, i.e. 220 GHz, the 
resolution of the 30-m telescope is about 11". The dump time was 2 s, and the sampling interval 
4". All the maps except one were scanned twice, once along the R.A. direction and once along the 
Dec. direction, to avoid scanning effects on the resulting images. The map of G19.61-0.23 was 
scanned only along the R.A. direction due to lack of time. The velocity resolution is ~0.1 km s ' 
(80 kHz), and the bandwidth ~1 10 km s ' (80 MHz). 

The data were reduced and analysed with the programs CLASS and GRAPHIC of the GILDAS 
software package developed by the IRAM and the Observatoire de Grenoble. Both the '^CO(2-l) 
and the C'^0(2-l) spectra have been smoothed to a resolution of 0.5 km s in order to improve 
the signal to noise ratio. After smoothing, the average rms level for a single spectrum is ~2.5 K, 
and the average rms of the outflow maps is ~0.3 K. 

3. Results 

3.1. Molecular outflows 

In Fig. [T]we present the '■^CO(2-l) (solid) and C'^0(2-l) (dashed) spectra for each observed 
source. These have been obtained by averaging the emission inside the half power contour of the 
C'**0 integrated maps (white solid contour in Fig.|2]l. A comparison between the '''CO(2-l) and 
the optically thinner C"*0(2-l) lines allows us to define the systemic velocity, Vlsr^ of the molec- 
ular clumps, identify the '^CO line wings and discriminate between self-absorption and multiple 
components in the '^CO(2-l) line profiles. 

All the '^CO(2-l) spectra display high-velocity wings, typical of outflowing motions. Outflow 
maps are shown in Fig.|2l where the blue- and the red-wing emission contours (solid and dashed, 
respectively) of the ^^CO(2-l) line are overlayed on the C'^0(2-l) integrated map (grey scale). 
H2O and CH3OH maser spots are also indicated on the maps as filled squares and triangles, respec- 
tively, and the positions of the UC Hii regions are marked as crosses. 

Out of the 1 1 sources observed, 9 show red- and/or blue-shifted lobes. The remaining two 
sources (G23.44 and G29.96) both show high velocity '^CO(2-l) wings and hence presumably 
have outflows also (in the case of G29.96, see discussion of Beuther et al. |2007bl ) but they are 
sufficiently confused with our angular resolution that we do not discuss them further here. Higher 
angular resolution maybe in a different tracer is needed to resolve any possible outflow in these 
regions. The red emission in G10.47 is dominated by a secondary nearby clump towards the south- 
west, marked as G10.46B in Fig.|2]because of its proximity to the UC Hii region bearing the same 
name (Cesaroni et al. |1994t see Appendix). As a consequence, a well-defined isolated red lobe 
does not appear on the map at the resolution of the observations. However, the elongated blue 
lobe, together with the wide wings of the '^CO(2-l) spectrum favour the presence of a molecular 
outflow in this region. 

3.2. Molecular clumps 

The good spatial overlap between the C"*0 bulk emission and the H2O and CH3OH maser spots 
(Fig. 12]) confirms that the C"^0 molecule is efficiently tracing the high-density clump hosting the 
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55 60 65 70 75 80 -15-10-5 5 10 50 55 60 65 70 




30 40 50 60 90 100 110 90 100 110 120 




5 10 15 20 25 30 -10 10 



Vlsr (km/s) 

Fig. 1. '^CO(2-l) (solid) and C'^0(2-l) (dashed) spectra of each source, obtained by integrating 
the emission inside the 50% contour level of the C'^0(2-l) maps. In each box the scale to the 
left corresponds to the '^CO(2-l) main beam brightness temperature, the one to the right that of 
C"^0(2-l). 



massive YSOs. The white solid line in each map represents the 50% contour level of the C'**0(2-l) 
integrated emission and gives a measure of the size of the molecular clump. In columns 2 and 
3 of Table |2] the deconvolved linear diameter of the C'^O clumps, D, and the optical depth of 
the C'^0(2-l) emission, tcisq, are presented. The latter has been estimated from the '^CO(2-l) 
to C'**0(2-l) intensity ratio measured at the systemic velocity, using the spectra in Fig. [T]and 
assuming the same excitation temperature, Tex, for both species. A relative abundance of 0.14 has 
been assumed between the C"*0 and '^CO isotopologues (Wilson & Rood |1994] l. 

From the emission integrated under the C^^O(2-l) line corrected for tcisq. the clump mass, 
A^ciump, has been determined assuming an excitation temperature, Tex, between 20 and 80 K, and a 
C'^O abundance dependent on the Galactocentric distance (Wilson & Rood 1 19941 ): 

[C'^O] 

1^ = + 37.1 (1) 
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Table 2. Deconvolved clump diameter and mass estimates 
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Source 


D 


TclSo" 






Tdust 






A^dyj'] 




(pc) 




(Mo) 




(Mq) 




(km s ) 


(Mq) 


G10.47 


1.25 


0.57 


1460-3200 


5680 


20' 


47900' 


0.57 


44 


G10.62 


0.98 


0.82 


1920-4110 


4770 


20' 


20200' 


0.61 


27 


G16.59 


Q.ll 


0.69 


390-870 


1210 


20' 


1130' 






G19.61 


1.19 


0.60 


2450-5160 


5320 


20' 


33000' 






G23.44 


1.79 


0.76 


2390-5120 


2850 


20' 


4160' 






G28.87 


0.96 


1.40 


1260-2790 


1350 


33^ 


3420^ 






G29.96 


1.34 


0.56 


2030-4460 


2030 


20' 


12900' 






G35.20 


0.32 


0.83 


110-250 


460 


20^ 


1500^ 


0.66 


14 


G43.89 


0.91 


0.85 


570-1200 


2230 


10" 


6610* 


0.86 


88 


G48.61 


0.83 


0.27 


1200-2620 


1280 


29^ 


3790' 






G75.78 


0.54 


0.41 


450-1000 


1390 






0.68 


23 



" Optical depth of the C"*0(2-l) emission, estimated from the '^^C0(2-l) to C'^0(2-l) line intensity ratio 
at the systemic velocity 

'' Clump mass estimated from the C'*0(2-l) emission using the area under the whole emission line, and the 
spatial area delimited by the half power contour. Tex assumed to be between 20 and 80 K 
Virial mass, computed using Eq. (6) of Fontani et al. C002b 

Clump mass derived from the (sub)mm dust emission measured by the cited authors 
Rotational velocity; projected values (i.e. lower limits) 
' Dynamical mass of the clump, estimated using Eq. l|3j 
References: ' Hill et al. ( 12005^ . ^ Faundez et al. ( |2004] |. ^ Mooney et al. ( fT995l i. Hatchell et al. ( |2000] |. ^ 
Mueller et al. ( l2002t 



We have adopted [CO]/[H2] - 10 Fractionation and selective photo-dissociation have been 
neglected. Table [3]presents the results from the gaussian fit to the C'^0(2-l) spectrum at the peak 
position of each integrated map (grey scale in Fig.|2]i. The line wings were ignored in the fit. The 
virial mass, Mvir, of each clump, has been determined from the FWHM thus estimated, and from 
the angular diameter of the clump, using Eq. (6) of Fontani et al. ( 120021 1. The values obtained for 
A^ciump and Mvii are listed in Table|2l It is evident, from a comparison between these two variables, 
that the virial masses are roughly of the order of the clump masses, indicating the molecular clumps 
are in virial equilibrium. 

For comparison with the C'^O and virial masses, we have taken from the literature observations 
mapping the mm or submm continuum flux and we have converted these (see col. 7 of Table|2]l into 
estimates of the clump mass inferred from the dust emission, Mdust- For this purpose, we used the 
equation: 

Mdust = (2) 

where 5 V is the measured flux, d the distance to the source (see Table[Tl), /Cy the frequency-dependent 
dust mass opacity coefficient, By(Td) is the Planck function for a blackbody of dust temperature Tj, 
and is the dust-to-gas ratio, which we assume equal to 0.01. We have assumed /Cy = Ko(v/vo'f, 
with Ko = 1 cm^ g"' at vq = 250 GHz (Ossenkopf & Henning |1994| l, and jS = 2, as in Beuther et al. 
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( |2002al l. Td varies from clump to clump according to the authors who measured the dust emission 
(see Table|2]l. We note that no (sub)mm observations towards G75.78 are available in the literature, 
and that the dust mass obtained for G48.61 has been derived from 350 fim data. The issue of clump 
masses will be further discussed in Sect. 5. 

Table 3. Results from gaussian fit to C'^0(2-l) peak spectra 



Source 


Tub 


FWHM 


Vlsr 


/ Tub dV 






(K) 


(kms-') 


(kms-') 


(Kkms-') 


GIO 


.47A 


13.1 (3.1) 


6.6 (0.3) 


66.4 (0.1) 


92.1 (3.8) 


GIO 


.47B 


7.4(1.0) 


4.9 (0.4) 


71.1 (0.2) 


38.6 (2.5) 


GIO 


.62 


31.3 (1.7) 


6.8(0.1) 


-3.2(0.1) 


227 (2) 


G16 


.59 


13.6(1.3) 


3.9 (0.1) 


59.0 (0.1) 


56.2(1.2) 


G19 


.61 


14.4(1.8) 


6.5 (0.2) 


41.5 (0.1) 


99.8 (2.5) 


G23 


.44 


14.0 (3.2) 


3.9 (0.3) 


100.8(0.1) 


58.0 (3.5) 


G28, 


.87 


15.6(1.6) 


3.7 (0.1) 


103.3 (0.1) 


60.9 (1.7) 


G29 


.96 


21.0(1.4) 


3.8(0.1) 


97.6 (0.1) 


84.9(1.4) 


G35 


.20 


10.1 (1.1) 


3.7 (0.1) 


34.9 (0.1) 


40.4(1.2) 


G43, 


.89 


8.0 (2.0) 


4.8 (0.4) 


53.9 (0.1) 


41.1 (2.4) 


G48, 


.61A 


12.2 (2.3) 


3.8 (0.2) 


18.6 (0.1) 


49.7 (2.4) 


G48, 


.61B 


3.4 (0.7) 


5.6 (0.6) 


15.8 (0.2) 


20.4(1.8) 


G75, 


.78 


10.7 (2.4) 


5.0 (0.3) 


-0.2 (0.1) 


56.8 (3.0) 



In some cases we detect velocity gradients in the C'^0(2-l) emission (Fig. [3). The interpre- 
tation of these is not straightforward, but the comparison between the direction of the velocity 
gradient and that of the main outflow axis can give some hints. If the velocity gradient is parallel 
to the outflow axis, it could be reflecting the molecular outflow akeady observed in the '^CO(2-l) 
emission. On the other hand, if this gradient is perpendicular to the outflow axis, as is the case 
of, e.g., G35.20, it could be interpreted as rotation, similar to what is found in G24.78+0.08 and 
G3 1.41 +0.31 (Beltran et al. 2004). 

Whichever the nature of the velocity gradients, it cannot be firmly confirmed with single-dish 
observations alone. Complementary interferometric observations are necessary to assess their na- 
ture. However, under the assumption that there are rotating structures, we have estimated a dynam- 
ical mass, Mdyn, for the five cases in which we detect a C^^O velocity gradient, from the equation: 

R 

Mdyn = -g- (3) 

where Vjot is the rotation velocity, R the radius of the supposed disk, and G the universal gravi- 
tational constant. This equation assumes equilibrium between centrifugal and gravitational forces, 
and yields only a lower limit to the actual dynamical mass because Vrot is not corrected for incli- 
nation. Both Vrot and Mdyn are presented in Table |2l Our values are small compared to the other 
clump mass estimates, and of the same order of magnitude as those found by Beltran et al. ( 120041 ) . 
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4. Derived outflow parameters 

4. 1 . Outflow velocity ranges 

The determination of outflow parameters from our '^CO spectra is non-trivial. The first compH- 
cation arises when trying to select the outflow velocity ranges: there are no universal methods to 
separate the outflow from the ambient gas. If one includes too low-velocities there might be contam- 
ination from the ambient gas and the mass of the outflow may be overestimated. On the contrary, if 
one prefers to be safe from including ambient gas emission and chooses only the highest velocities, 
an important part of the outflow mass may be missed. 

The strategy we have used to select the blue and red velocity ranges of the outflows is the 
following. From the '^CO(2-l) and C'^0(2-l) spectra in Fig.[Tl we have defined the low-velocity 
limits where the C'^0(2-l) line reaches a level of 2cr or, in some cases, the beginning of weak 
non-gaussian wings. The high-velocity limits have been chosen where the '■'CO(2-l) falls below 
2cr. With the extra aid of the channel maps, these limits were either conserved or slightly modified 
depending on the blue- and red-shifted emission spatial distribution observed. The resulting blue 
and red velocity ranges, AVbiue and AVied, are listed in Table H] These ranges have been used both 
to make the outflow maps in Fig.|2]and to calculate the outflow parameters. For G10.47, no AVred is 
shown because the red-shifted outflow emission is contaminated by a second clump located ~30" 
towards the south-west. Similarly, no AVbiue is given for G28.87 because the blue emission does 
not show a well-defined structure in the map. 

Table 4. Velocity range of the '■'CO(2-l) wings used for the calculations and lobe sizes 



Source 










Sizebiue 


Size^d 




(km s-') 


(km s-') 


(km s-') 


(km s-') 


(pc) 


(pc) 


G10.47 


[56.3, 60.8] 




10.1 




24 




G10.62 


[-11.6, -8.6] 


[3.9, 7.9] 


8.4 


11.1 


0.9 


0.7 


G16.59 


[55.0, 56.5] 


[62.0, 64.0] 


4 


5 


0.5 


0.5 


G19.61 


[29.1, 35.1] 


[48.6, 50.6] 


12.4 


9.1 


0.3 


1.2 


G28.87 




[108.1, 113.6] 




10.3 




0.8 


G35.20 


[25.3, 30.3] 


[37.8, 43.8] 


9.6 


8.9 


14 


0.8 


G43.89 


[41.7, 49.7] 


[57.2,61.2] 


12.2 


7.3 


0.3 


0.4 


G48.61 


[9.7, 11.2] 


[21.7, 25.2] 


8.9 


6.6 


0.8 


0.9 


G75.78 


[-9.1, -4.1 


[4.4, 10.9] 


8.9 


11.1 


0.5 


0.4 



4.2. Optical depth 

Thanks to the high signal-to-noise level of the observations in the '■'CO(2-l) wings and even 
in most of the C'**0(2-l) line wings, it has been possible to estimate the optical depth of the 
'^CO(2-l) outflow emission (which we denote as Tbiue and Tred for the blue and red emission, 
respectively) as done for Tc'so (see Sect. 3.2). One should realise that the S/N ratio in the C'*'0(2- 1 ) 
wings is sufficiently high only for the lowest-velocity channel of each outflow velocity range, and 
therefore Tbiue and Tred have been computed only for such channels. Consequently, the values we 
obtain should be taken as upper limits to the mean optical depth throughout the '^CO(2-l) wings. 
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and the same applies to the derived outflow parameters (mass, momentum, etc). In the case of 
G28.87 and G43.89, the signal in the C''^0(2-l) line wings was not high enough to provide an 
optical depth estimate. In Table |5] we do not indicate Tbiue and Tied, but the mean value, Tm - 
(^biue + T"red)/2. We notc, howcvcr, that even though we will refer to in the text from now on, it 
is actually Tuue and rred that have been used in the calculations. 



4.3. Outflow energetics and kinematics 

For each outflow, the mass, momentum and energy (M, p and E, respectively) have been derived 
from the '^€0(2-1) emission under each line wing, AVbiue or AV^ed, as follows: 



M ^YjiHch (4) 

ch=l 

n 

P = ^ MchlVch - Vsysl (5) 

ch=l 
1 " 

ch=l 

where the summation is for all the velocity channels within AVbiue or AVred, Vd, is the velocity of 
the corresponding channel, Vsys the systemic velocity of the clump (from Table |3j, and nid, is the 
mass measured in each channel within the area defined by the 5cr contour level of the corresponding 
outflow lobe (i.e. the most external outflow contours in Fig.|2]l. The sum of the contributions from 
each pixel and velocity channel builds up the total outflow mass, momentum and energy. The 
abundance ratio ['^CO]/[H2], as in the case of C'^O, depends on the Galactocentric distance of 
each source according to the following relation (Wilson & Rood |1994] l: 

[t3^ = 7.5^,3, + 7.6 (7) 

The excitation temperature assumed is Tex = 10 K (the numbers decrease by ~ 10% if Tex = 20 K 
is used instead of 10 K, and coincide for Tex - 30 K). The total outflow mass, momentum and 
energy obtained are presented in Table|5j both for the case r « 1 and corrected for Tm. On average, 
Tjn - 2, which is not negligible and causes M, p and E to increase by a factor ~2.5. 

A kinematic timescale, fjdn, for each flow lobe has been determined from the equation fjdn - 
f/Vom, where r is the projected size of the outflow lobe (see Table|4|, and Vout is the difference in 
absolute value between the flow velocity and the systemic velocity (Vmax in Table|4]l. The numbers 
thus obtained for the blue and red lobes have been averaged to produce the kinematic timescales 
given in Table |6] 

From fkin, the outflow parameters M, p and Lmec, have been derived for the blue and red lobes 
separately, and then added together to obtain the total values presented in Table |6l 



4.4. Outflow inclination 

We have no means to determine from our data the inclination, of each outflow. Here we assume 
i as the angle between the outflow axis and the line of sight. Therefore the parameters we have 
derived so far (Tables|5]and|6l) are subject to an uncertainty. As an indication, in Table|2]we show the 
range of values for the outflow parameters after correcting for ; of between 30°and 60°. Quantities 
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Table 5. Total outflow energetics. 







M 




P 












(Mo) 




> 


f 1 fl** ( 


erg) 






T <C 1 




T <K 1 




T <C 1 




n 1 n Alt 




150 


570 


1 1 1 n 


zLi 


o. 1 


JV 




1 Q 


yo 


zUU 


D / U 


1 s 1 n 


J.U 


1 9 


ni A -^Q 
oio. jy 


1 1 


24 


59 




1 1 n 
ziu 


U.J 1 


U. / J 


G19.61 


3.6 


200 


740 


1740 


6320 


15 


55 


G28.87* 




70 




470 




3.3 




G35.20 


1.6 


110 


230 


610 


1220 


3.5 


6.9 


G43.89 




62 




380 




2.5 




G48.61 


0.9 


590 


900 


2550 


3870 


11.8 


18 


075.78 


0.9 


140 


220 


920 


1370 


6.0 


9.0 



^ Mean t of the ''C0(2-l) line over the blue and red wings 

* Only the blue lobe parameters, for G10.47, and the red lobe parameters, for G28.87, were determined, 
due to contamination from nearby unrelated clumps 

Table 6. Total outflow timescales and kinematics. 



Source 




4in 




M 




P 










(lO-^yr) 


(10-3 


Moyr-') 


(10-3 Mq 


km S-' yr-') 


(L©) 










T « 1 


Tm 


r « 1 




T « 1 


Tm 


G10.47* 


3.6 


23 


0.66 


2.4 


4.7 


18 


2.9 


11 


G10.62 


1.9 


8.4 


1.1 


2.5 


8.4 


20 


5.5 


13 


G16.59 


2.1 


11 


0.23 


0.56 


0.82 


2.0 


0.26 


0.62 


G19.61 


3.6 


8.0 


4.4 


17 


38 


150 


28 


110 


G28.87* 




8.1 


0.85 




5.8 




3.4 




G35.20 


1.6 


11 


1.0 


2.0 


5.2 


10 


2.5 


5.0 


G43.89 




4.0 


1.9 




12 




7.0 




G48.61 


0.9 


11 


6.1 


9.3 


26 


39 


10 


15 


G75.78 


0.9 


4.5 


3.5 


5.2 


22.3 


33 


13 


18 



"I" Mean t of the '3CO(2-l) line over the blue and red wings 

* Only the blue lobe parameters, for G10.47, and the red lobe parameters, for G28.87, were determined, 
due to contamination from nearby secondary clumps 

that depend strongly on the flow velocity, i.e. the kinematic energy and the mechanical luminosity, 
are most affected by the inclination, especially at high values of /. The numbers variate from a 
factor -0.6 for M at 30°to a factor ~7 for L^ec at 60°. 

5. Discussion 

5.1. '^CO versus '^CO 

The '^CO molecule is about 77 times more abundant than '^CO in the interstellar medium (Wilson 
& Rood 1 19941) . and consequently optically thicker In order to understand the effects of using one 
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Table 7. Total outflow parameters after correction for an inclination angle in the range 30° - 60° 
(Tm taken into account) 



Source 


D 
t 


E 


M 


t 






(Mo km s"') 


W'^ erg 


(lO^'Mo yr"') 


(10"' Mo km s"' yr"') 


(Lo) 


G10.47 


4770-8260 


40-120 


\A-A.1 


12-61 


8.3-75 


G10.62 


1740-3010 


15-46 


1.4-4.3 


13-68 


10-92 


G16.59 


240-410 


1.0-3.0 


0.32-0.97 


1.3-6.8 


0.48-4.3 


G19.61 


7300-12600 


73-220 


9.8-29 


97-500 


82-740 


G28.87'" 


540-940 


4.3-13 


0.49-1.5 


3.9-20 


2.6-24 


G35.20 


1410-2400 


9.2-28 


1.1-3.4 


7.0-36 


3.8-35 


G43.89t 


440-770 


3.3-10 


1.1-3.3 


8.0-42 


5.3-48 


G48.61 


4470-7740 


24-71 


5.3-16 


26-140 


12-100 


G75.78 


1590-2750 


12-36 


3.0-9.0 


22-110 


14-130 



+ G28.87 and G43.89 have no t„, estimate 

tracer or another on the derivation of outflow parameters, we have made the following analysis: 
we have compared our outflow parameters, obtained from our '^CO(2-l) observations, to those 
derived by other authors from '^CO for the sources G16.59-0.06 and G75.78+0.34. The former 
belongs to the sample surveyed by Beuther et al. ( |2002bl l and has been studied in detail by the same 
authors more recently (Beuther et al. I2006I I. both with IRAM 30-m and with the Sub Milhmeter 
Array (SMA); the latter has been mapped at high angular resolution by Shepherd et al. ( 119971 see 
Appendix). Interferometric observations were also performed towards G16.59 in the '^CO(l-O) 
transition by Furuya et al. ( I2008I I, revealing just one of the two outflows resolved by Beuther et 
al. ( I2OO6I 1 with a combination of SMA and IRAM 30-m images, and so we do not use their results 
in this comparison. To avoid effects due to the uncertainties in the kinematic timescales of the 
outflows we use only M, p and E in the comparison, which are time-independent. 

We note that Beuther et al. ( |2006l l chose the same low-velocity limits as us in AVbiue and AVied, 
but used a '^CO abundance of 1.1x10"^ relative to H2, which in the case of G16.59 is smaller 
than the one we have used by a factor 2.3. Therefore, just for the purpose of this section, we 
have modified our G16.59 outflow parameters accordingly. On the other hand. Shepherd et al. 
( I1997I I chose low-velocity limits higher than ours, missing the lowest-velocity channels that we 
have included in our calculations. We have therefore re-computed our outflow parameters using 
the same limits as Shepherd et al. ( 119971 ) to allow for a coherent comparison. The numbers are 
presented in TableHl together with the tracers used by the cited authors. 

From this table it is clear that, when the AV chosen are the same, all our parameters are compa- 
rable to those derived by Beuther et al. ( I2006I I and Shepherd et al. ( 119971 ). Our mass estimates are 
somewhat higher, even when t « 1 is assumed. This suggests that the optical depth of the '^CO 
line wings may have been underestimated by the cited authors. It is also true that our r,n have likely 
been overestimated for the reasons described in Sect. 4.2. That is why we also present our results 
for the optically thin case in Tables |5] and |6] constraining our outflow parameter estimates to fall 
between the two limits. 

On the other hand, our energy estimates are smaller on average than those derived by Beuther 
et al. ( I2OO6I 1 and Shepherd et al. dl 9971 1. This is due to the more extended wings of the optically 
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Table 8. Comparison of outflow mass, momentum and kinetic energy between this work and other 
studies 



Source 


Tracer 


Mother 






Pother 


Pour 




^other 










(Mq) 


(Mo) 




(Mokms-') 


(Mq km 




(10"' erg) 


(10^^ er 


g) 








T 1 






r « 1 


Tra 




T « 1 




G16.59^ 


'-C0(2-l) 


33 


55 


130 


268 


190 


480 


2.5 


0.71 


1.7 


G75.78* 


'-CO(l-O) 


54 


66 


100 


667 


520 


800 


8.3 


4.1 


6.3 



^ Another, Pother and Mother froui siugle-dish observations by Beuther et al. ( I2002bll2006b . corrected for revised 
distance 

'" ^oiher> Pother fud Mother from Interferometric observations by Shepherd et al. J1997| l, corrected for revised 
distance; our outflow parameters re-calculated using the same low-velocity limits as Shepherd et al. ( 1 19971 1 



thicker '^CO Unes. Unlike the mass, the kinetic energy estimate is most sensitive to the higher 
velocity channels (see discussion by Masson & Chernin |1994| l. and so it is to be expected that a 
line with broader wings will give rise to higher kinetic energies. Moreover, the effect of the optical 
depth is much smaller in this case because it is negligible in the high-velocity channels. 

A compromise between the mass and the energy is the momentum of the flows, which does 
not depend so much on the extremes of the velocity ranges chosen (again, see Masson & Chernin 
1 19941 1. Indeed, we see from Table |8] that p is the parameter that shows the smallest discrepancy 
between different studies. 

We conclude from this comparative analysis that the use of the '^CO isotopologue will gener- 
ally produce more reliable masses, similar momentum and underestimated kinetic energies relative 
to what is obtained using '^CO. This effect is evident in Fig.|4](see Sect. 5.2) 

Note that the loss of the highest velocity channels when using '^CO instead of '^CO also leads 
to higher kinematic timescales. Indeed, the values of fun obtained in the present work are on average 
~1.5 times higher than those derived by Beuther et al. ( |2002bl ), a result which can be explained by 
the fact that our mean outflow Vmax is a factor ~1.5 lower. 

5.2. Outflow parameters as a function of luminosity 

In this section, we discuss our results and compare them to those obtained by Beuther et al. (I2002bl) . 
who carried out a similar survey in less luminous sources. 

In Fig. m the time-dependent outflow parameters, i.e. mass loss rate, M, mechanical force, p 
or Font, and mechanical luminosity, Lmec, are plotted against bolometric luminosity, Lboi (Table[T]i. 
The quantities obtained by Beuther et al. (I2002bl l. for their sources with no distance ambiguity, are 
represented as open circles in Fig.|4] and our values as filled symbols: circles for the sources with a 
Tm estimate, and triangles for G28.87 and G43.89, for which no optical depth has been determined 
in the '■'CO(2-l) line wings. The latter represent lower limits to the outflow parameter values. 
The same symbology applies to Figs. |5] |2]and|8] The error bars correspond to a correction for an 
inclination varying from 30° to 60° (see Table IT). It is immediately apparent from these plots that 
our sources have higher bolometric luminosities than those studied by Beuther et al. by a factor 
~10, and thus complement their study. 
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The results by Beuther et al. ( |2002bl l show continuity in the correlation between mechanical 
force, f out, and bolometric luminosity, Lboi, from low-mass to high-mass outflows, with higher 
mechanical forces corresponding to higher luminosities. In this work we find that this trend extends 
up to luminosities higher than 10^ L0, indicating that the outflow formation mechanism is present 
up to luminosities of the order of 10* Lq, and that more luminous sources drive more powerful 
outflows. Similar correlations are found for the mass loss rate, c/Mout/t/f, and for the mechanical 
luminosity, Lmec, of the outflows. For the latter, there seems to be a saturation towards the highest 
luminosity sources, but this may just be a consequence of using '^CO instead of the '^CO molecule 
used by Beuther et al. (|2002bl as explained in Sect. 5.1). The weakening of the correlation between 
dMout/dt and Lboi towards the highest luminosities claimed by Beuther et al. (I2002bl l is not evident 
in our plot: higher luminosities yield higher mass loss rates, with no apparent saturation in the 
highest luminosity range. 

The behaviour of the data plotted in Fig.|4]is consistent with the results obtained by Wu et al. 
(120041 Figs. 6 and 7: 120051 Fig. 13), who compiled all the outflow studies available before February 
2003 and performed a statistical analysis of the outflow parameters for a sample of sources spanning 
a wide range of bolometric luminosities (10 '-10* L0). Despite the heterogeneity of the sample, 
their best fits for Mout, ^out and Lmec (vs Lboi) fit nicely the points in Fig. |4] (dashed lines). As 
pointed out by the cited authors, the fact that these correlations hold for such a broad range of 
luminosities suggests that the luminosity of the powering source determines the outflow energetics, 
and that the driving mechanisms are similar for all luminosities. 

It is unclear from our observations which is the multiplicity of the outflows. The increasing Font 
with luminosity could be interpreted both as a single outflow powered by a more luminous source, 
or by more than one outflow coming from less luminous sources, each of them contributing to the 
overall outflow energetics (see Appendix for some particular cases). Such outflow multiplicity is 
generally found to be around 3 or 4 when interferometric observations are available (e.g. Beuther 
et al. l2006l l. If we assume that the highest luminosity clumps contain a few molecular outflows, the 
only effect on the plot would be to split each point into several points shifted by a small factor (3 
or 4) towards lower Fout and Lboi, so the overall appearance of the plots in Fig.|4]would remain the 
same. 

For completeness, we also present plots of the time-independent outflow parameters as a func- 
tion of Lboi (Fig-El), which show an analogous behaviour to that of Fig.H] We feel that, based on 
the arguments described, the general behaviour of the outflow properties plotted in Figs.|4]and|5]is 
real and not affected by possible errors in the calculation of the outflow parameters. This belief is 
reinforced by the continuity displayed by our data with respect to that of Beuther et al. (I2002bl l. 

In Fig. |6] we plot the Lyman photon rate of the associated UC Hii region(s) (see Table [1]), 
A^Ly, against bolometric luminosity, Lboi. We have computed A^Ly for the sources from Beuther et 
al. ( |2002bl ) which were detected in the VLA 3.6 cm survey by Sridharan et al. ( !2002| l, under the 
assumption that all the emission comes indeed from an Hn region and that it is optically thin. The 
dashed line depicts the expected curve for an individual source on the ZAMS (Panagia |1973l ). The 
fact that most of the sources lie below the theoretical curve suggests that the observed regions 
contain a whole cluster of YSOs of different masses. All of them contribute to the total Lboi, but 
only the most massive stellar objects, which power Hii regions, contribute to N^y. G43.89 lies just 
above the theoretical ZAMS line, but as indicated by the arrow, an increase in the distance to this 
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source would shift the point below the line. The same applies to the two other sources from Beuther 
et al. (I2002al) sample which are also above the ZAMS line. It therefore becomes evident that an 
accurate distance determination is necessary in order to interpret the observations correctly. 

In an attempt to assess our data as a function of a quantity which does not represent a whole 
cluster, such as Lboh but the individual massive stellar objects, in Fig. |7]we plot Fout against the 
Lyman photon rate. There seems to be a correlation between Font and A^Ly which holds within a 
range spanning from B2 to 06 stars, assuming the UC Hii regions are powered by a single stellar 
object. 

Despite being less affected by source multiplicity, this plot displays a similar correlation as 
FoutiLhoi) in Fig.m and indicates that more massive objects power more energetic outflows. Even if 
there are multiple UC Hn regions contributing to the total ionising photon rate, the good correlation 
between Fout and A^Ly shown by this plot would suggest that each outflow can be associated with 
an individual massive stellar object. 

5.3. Outflow parameters versus clump mass 

Another way of inspecting our data is to plot some outflow parameters as a function of clump 
mass. Figure [8] shows outflow mass. Mom, and mechanical force, Fo^t, against clump mass. On the 
left panels, our clump mass is the one obtained from the C'^0(2-l) emission (see Table |2]l, while 
masses of Beuther et al. (I2002bl) were determined from the 1.2 mm continuum emission of the 
clumps (Beuther et al. |2002al ). Therefore, for consistency in the comparison, we show on the right 
panels the same plots, with the only difference that our clump masses here correspond to Mdust in 
Table |2] (see Sect. 3.2). In order to homogenize the two samples, we have re-computed Beuther's 
clump masses from their measured 1 .2 mm fluxes (Beuther et al. I2002al ) using Eq. and we 
find them to be smaller by a factor 4.6 with respect to those reported in the original paper This 
difference is due both to our lower assumed grain opacity with respect to the one used by Beuther 
et al., and to an error in the grain emissivity approximation by the same authors (see their Erratum, 
published in 2005). Once more, it is evident that our data complement the data from Beuther et al. 
( |2002bl l, adding higher mass molecular clumps. 

Two main ideas can be extracted from an overall look at the left and right panels in Figure 
[8] Firstly, our dust clump masses are larger than the corresponding C'^O masses by an average 
factor of ~5. This difference might be explained by the fact that C'^O and the (sub)mm continuum 
are tracing different parts of the clump. Indeed, different sizes of the clump are measured in each 
tracer: the angular FWHM measured in the (sub)mm surveys is larger than in our C'**0(2-l) maps 
by a factor of ~ 2.5. This may account for great part of the difference between the two mass 
estimates. In addition, Hofner et al. (120001 ) concluded from their survey that the masses derived 
from (sub)mm dust emission tend to be systematically higher than masses derived from C'^0(2-l) 
and C'^0(2- 1) lines by a factor 2. They point out that several sources of uncertainty may contribute 
to this discrepancy, such as C^^O abundance, optical depth estimates, and the dust grain emissivity 
adopted. 

The second feature to note in Fig. [8] is the better continuity from Beuther's data towards our 
higher mass data when Mdust is used (right panels). This continuity is not so clear when we plot our 
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data using the C'**0(2-l) clump mass instead (left panels), especially for the Mout plot. This is not 
surprising, since we are comparing parameters derived from different tracers. 

We find that when Mjust is considered, one sees a tight correlation between this and the mass 
of the outflow, Mout- The approximate relation Mout - ^-^^dusv whi''^ is flie best fit to the data of 
Beuther et al. (I2002bl l. is also a fair fit to our more massive molecular clumps (see Fig. [8]l. 

In Fig. |9] we present a plot of Mout/Mjust against Lboi/Mdust, which are distance-independent 
quantities. Ignoring G48.61, which is the only source for which the dust mass has been estimated 
from 350 yum data (see Sect. 3.2), two main interesting features characterise this plot. Firstly, our 
points tend to cluster to the left of the points representing Buether et al. This indicates that the most 
massive clumps tend to be "under-luminous". Secondly, and most importantly, the points span only 
a factor 10 in Mout/Mdust, independently of their luminosity. This would suggest, consistently with 
Fig. [8] that the outflow mass is more strictly related to the mass of the clump than to the luminosity. 
In other words, the mass entrainment of the outflow is proportional to the mass of the whole clump. 

5.4. Accretion rates 

It is not straightforward to derive mass accretion rates from outflow mass loss rates. Nevertheless, 
by making a few assumptions it is possible to provide a rough estimate of the average accretion rate, 
Maccr, for our sample. Following the reasoning described in Beuther et al. (I2002bb . we consider the 
molecular outflows to be momentum driven. Therefore: 

Pout - Pjet = MjetVjetfkin (8) 

Here we have assumed that the kinematic timescale of the jet is the same as that of the molecular 
outflow (listed in Table |6]l, which in our case is of order 10^ yr The average pout, (from Table H) 
corrected for a mean outflow inclination of 45°, is ~4000 Mq km s ' . Adopting a typical jet veloc- 
ity of about 700 km s ' (e.g. Martf et al. 1 19981 1. the mass loss rate of the jet, Mjet, is estimated to be 
on the order of 6x10"^ M0 yr ' for sources with Lboi in the range 10^-10^ Lq. If we assume that 
Mjet is about one third of the accretion rate (Tomisaka |1998l ). Maccr amounts to ~2xl0 M0 yr 
which is consistent with other values previously reported for high-mass SFRs (e.g. Beuther et al. 
I2002bl Zhang et al. l2005T l. 

5.5. Velocity gradients 

For one of our sources displaying a velocity gradient in the C"*0(2-l) emission, and assuming 
such gradient is due to rotational motions, it is possible to find out whether the specific angular 
momentum, I, is conserved at dififerent scales. This source is G10.62, for which other studies have 
revealed a rotating structure at a sub-parsec scale. G35.20 has also been reported to contain a 
rotating envelope by Little et al. ( |19851 l. However, the size of such rotating structure is roughly the 
same as the one found in this work (-0.2 pc), and so we are not able to assess the conservation of I 
at smaller scales. 

Using Vrot in Table |2] and a radius of 0.6 pc, we obtain / ^ 0.2 km s"' pc for G10.62. On a 
smaUer scale (~0.1 pc), we derive I ^ 0.3 km s ' pc from the C'**0 interferometric observations 
carried out by Ho et al. (|1994| l. which is of the same order as the value obtained from our data. One 
can compare this value of / with the initial value of the parental clump. Assuming that the initial 
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clump angular momentum is due to the Galactic differential rotation, we find an angular velocity 
CO ~10"'^ s"' (Clemens |1985l l. For a 1000 M0 clump formed from a ~10 pc radius sphere with a 
density of 1 cm"^, typical of the interstellar medium, the initial I is hence ~3 km s ' pc, an order 
of magnitude higher than the derived (sub)parsec-scale value. 

To sum up, although this analysis is based only on one source, it seems that the specific angular 
momentum is roughly conserved at scales smaller than about 1 pc, but it is not conserved on larger 
scales. A similar result was found by Ohashi et al. ( 119991 see their Fig. 2) for a sample of low- 
mass protostars. In their case, the radius inside which / is conserved is smaller by about an order 
of magnitude, which indicates that infall takes place within a smaller radius than in the case of 
high-mass clumps. 

6. Summary 

We have mapped 1 1 high-mass SFRs in '^CO(2-l) and C'^0(2-l) with the IRAM-30m telescope 
(Spain), with the aim to search for molecular outflows and characterise them and their associated 
molecular clumps. The sample is composed by high luminosity (up to ~10* Lq) massive young 
clumps in their earUest evolutionary phases. 

Our main conclusions can be summarised as follows: 

1. The whole sample shows high-velocity wings in the '^CO(2-l) spectra, indicative of outflow- 
ing motions. We have obtained outflow maps in 9 of our 1 1 sources, which display well-defined 
blue and/or red lobes. This result indicates that outflows are as common in high-mass SFRs as in 
low-mass SFRs. Furthermore, it suggests that massive star formation may commonly proceed 
by accretion, in a similar way as low-mass star formation. 

2. From the C'*0(2-l) emission, the clump masses have been derived, taking into account the 
optical depth. The masses obtained are high and comparable to the virial masses, but have 
been found to be systematically smaller than the coiTesponding masses computed from the mm 
continuum emision by a factor ~5. This is largely due to the different sensitivities in the two 
tracers. 

3. The outflow parameters have been determined from the '^CO(2-l) emission in the high- 
velocity wings, by using standard methods. The values obtained are coiTected for optical depth, 
which we have been able to estimate from the '^CO(2-l) to C'^0(2-l) ratio in the line wings. 
Our results complement those by Beuther et al. (I2002bl) . adding more luminous sources by 
an order of magnitude on average. A comparison between our results and those of Beuther et 
al. (I2002bl) reveals that there is continuity in the behaviour of outflow-related quantities as we 
move to the higher luminosity sources in our sample: higher outflow-related values coiTespond 
to higher luminosity sources. This trend also agrees quantitatively with the correlations derived 
by Wu et al. ( 120041 12005I I for a large sample of outflow sources spanning a wide range of 
luminosities (O.l-lO*" L©). 

4. The mass of the outflow and the clump mass are tightly correlated, and the approximate relation 
^out - 0-3A^J^iump' obtained by Beuther et al. ( I2002bl ) fits adequately the observations over three 
orders of magnitude. 

5. We find a coiTelation between the outflow mechanical force, Font, and the rate of ionising pho- 
tons, A^Ly, of the UC Hii regions associated with the molecular clumps. High-angular resolution 
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images are needed to resolve the outflows and associate them with individual sources within 
the clump. 

6. Considering the molecular outflows in our sample to be momentum driven, we have estimated 
a mean accretion rate of ~ 2 x 10 Mq yr which is on the order of what is typically found 
in high-mass SFRs (e.g. Beuther et al. |2002bl l. 

7. We find C"*0 velocity gradients in 5 sources, in some cases they are compatible with rotat- 
ing envelopes akeady detected with high-density tracers and at high-spatial resolution by other 
authors. The detection of velocity gradients on such large-scales with the C'^O molecule en- 
courages further observations towards these five sources at higher- angular resolution and with 
high density tracers to understand their true nature. 
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Appendix A: Comments on individual sources 

A.1. G1 0.47+0.03 

With Lboi - 139 X 10* L0, this is the most luminous source in our sample. This region con- 
tains 4 UC Hii regions, three of which are concentrated inside an area of 2" across inside a hot 
molecular core (Cesaroni et al l 19981 1. coinciding with H2O masers (Hofner & Churchwell, 1996). 
High-angular resolution C'^O maps obtained by Gibb et al. ( |2004| l reveal a velocity gradient in the 
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SW-NE direction, which the authors interpret as a group of smaller individual cores with slightly 
different velocities. They conclude that there is weak evidence for either outflow or rotation in this 
region. On the other hand, high-resolution observations with PdBI by Olmi et al. (|1996l l show an 
east to west CH3CN velocity gradient, which they believe to be a rotating disk-like structure, and a 
south to north '^CO(l-O) velocity gradient which, added to the wide wings displayed by the line, 
could represent a molecular outflow. 

On a larger scale, Klaassen & Wilson (120081 ) report an HCO^ outflow in the direction NE-SW, 
with a linear scale of ~4 pc at 10.8 kpc. However, they stress that higher-resolution imaging is 
needed to confirm or reject this interpretation. In fact, from a close inspection of our maps we 
find that what they interpret as a molecular outflow is most likely a pair of cores with different 
systemic velocities (see Fig. IA.lt . The contamination by the second core is actually the reason why 
we cannot obtain a clear red outflow lobe map from our single-dish observations. 

As for the C'^0(2-l) first moment, the velocity gradient seems the least reliable of the five 
shown in Figure[3] Certainly the use of higher density tracers would be more adequate to assess the 
presence and nature of this velocity gradient. 

A.2. G10.62-0.38 

This source has been extensively studied in the past. It hosts an UC Hn region embedded in a hot 
molecular core (Keto et al. 119881 Plume et al. |1992| l. Its high Lyman photon rate (see Table ^ 
implies an 06 spectral type for a ZAMS star, if the Hn region were powered by a single stellar 
object. However, VLA radio continuum maps performed by SoUins et al. (2005) reveal a clumpy 
structure of the ionised gas, suggesting the presence of multiple stars. 

Several authors have concluded from NH3 observations that the HMC is rotating and collapsing 
inward onto and through the UC Hii region, and that rotation is seen at scales from 1 pc down 
to ~0.05 pc (e.g. Sollins et al. 1 20051 and references therein). The direction of the corresponding 
velocity gradient found by these authors goes from southeast (blue) to northwest (red) at small 
scales (2"). Our results show roughly an east to west velocity gradient in the C"*0(2-l) line (Figure 
[3]l, which despite having a slightly different orientation, could represent the large-scale continuation 
of the rotating core seen at high resolution by the above mentioned authors. Several H2O maser 
spots are aligned along the direction of the C'^0(2-l) velocity gradient. 

Our '^CO(2-l) outflow map displays well defined lobes and suggests an outflow axis along a 
northeast-soutwest direction, perpendicular to the rotating core seen from high resolution observa- 
tions (e.g. Sollins et al. 120051 Keto et al. I1988I I. The outflow axis coincides with the direction of 
the ionised jet discovered by Keto & Wood (I2006I I. 

A.3. G16.59-0.06 

Our outflow maps (Fig.|2]) suggest the presence of a bipolar outflow oriented in the NE-SW direc- 
tion, with a compact blue lobe, and a more extended, less collimated red lobe. 

Detailed studies of this region have been performed by several authors. In particular, Furuya 
et al. ( I2OO8I 1 obtained high-angular resolution '^CO(l-O) maps with the PdBI, which reveal a SE- 
NW molecular outflow. Beuther et al. (120061) . who mapped G16.59 in the '^CO(2-l) fine both with 
IRAM 30-m and with the SMA, resolve two perpendicular molecular outflows by combining their 
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single-dish and interferometric images: one flow is oriented along the same direction reported by 
Furuya, and the other along the NE-SW direction. Our observations are therefore consistent with 
the results obtained by Beuther et al. ( |2006l l. 

A.4. G35.20-0.74 

The presence of a bipolar molecular outflow in this clump is strongly suggested by the elongated 
structure of the blue- and red-shifted '^CO(2-l) emission (Fig. |2]i. Although the blue and red 
outflow lobes overlap each other, an outflow axis can clearly be defined along the SW-NE direction, 
and the overlap could be due to an outflow inclination almost perpendicular to the line of sight. 
This bipolar outflow has been previously observed in CO and SiO by Gibb et al. ( I2003I I. The 
same authors also present interferometric radio continuum maps showing a radio-jet oriented in 
the north-south direction, which is also detected in the mid- and near-IR by De Buizer (I2006I I 
and Fuller et al. ( 120011 ), respectively. Gibb et al. (120031 ) argue that the region contains up to four 
outflows, and consider this interpretation more plausible than precession of the jet to explain such 
different jet/outflow orientations. 

The C'^0(2-l) half power contour (white solid line in Figure|2]) has an elongated shape which 
is approximately perpendicular to the ^^CO outflow axis. In addition, a NW-SE velocity gradient 
is detected within this contour (Fig. [3]), which favours the interpretation of a flattened structure ro- 
tating perpendicularly to the outflow axis. Evidence of a large-scale disk-like structure was already 
presented by Little et al. ( 119851 ), who performed NH3 observations at various angular resolutions 
and detected a flattened molecular structure with the same orientation and size as our C'**0(2-l) 
velocity gradient. 

Our observations towards G35.20 are therefore consistent with results found by other authors, 
and are compatible with a rotating structure seen almost edge-on and perpendicular to a bipolar 
molecular outflow lying close to the plane of the sky. 

A.5. G43.89-0.78 

To our knowledge, no detailed studies on the kinematics of this high-mass SFR have been made in 
the past. From the wide wings in the '^CO(2-l) line (Fig.lT]) and the spatial distribution of the blue 
and red lobes in Fig. |2l we conclude that a high-mass molecular outflow is present in this clump 
(see also Tables |5] and |6|. 

Both the C'^0(2-l) peak position and the centre of the molecular outflow are offset with re- 
spect to the position of the associated UC Hn region (Wood & Churchwell 1 1 9891 ) by about 10", 
or ~ 0.2 pc at a distance of 4.2 kpc. This suggests that the outflow might not be driven by the 
massive YSO within the UC Hn region but by a different, presumably high-mass, young (proto)star 
embedded in the C'^O core. 

From the C'^0(2-l) first moment, a velocity gradient can be seen which is not as clear as in 
G10.62 or G35.20, but for which a southwest (blue) to northeast (red) direction can be assigned. 
This velocity gradient is neither aligned with nor perpendicular to the outflow axis (Fig. [3]), so 
its interpretation is unclear. Further observations at high-angular resolution would be helpful to 
understand the nature of the '^CO and C'^O emission. 
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A.6. G75.78+0.34 

Shepherd et al. ( 11 9971 1 conducted '^CO and SiO high-resolution observations towards this source 
with BIMA, and the resulting maps resolve three CO bipolar molecular outflows and one possible 
SiO outflow within an area of approximately 1' x 1' (roughly the area of the corresponding map 
in Figure |2]l. The axis of the reported CO outflows have a position angle (P. A.) ranging between 
north-south and east- west. Thus, our outflow map (Fig.|2]i, which shows no clear spatial separation 
between red and blue lobes, is probably englobing all three outflows detected by Shepherd et al. 
This is yet one more example in which the usefulness of complementary interferometric observa- 
tions is evident. 

More enigmatic is the presence of a C'**0(2-l) velocity gradient in this clump (Figure |3]l. Its 
west-east orientation is reversed with respect to two of the CO outflows discovered by Shepherd et 
al. (1997), while it is perpendicular to the third north-south outflow reported by the same authors. 
This would rule out the possibility that the velocity gradient is tracing an outflow, leaving us with 
the options of either rotation of the whole C"*0 clump or multiple components within it (or both). 
The former would imply a dynamical mass of 23 M0, too small to account for all the sources that, 
according to Shepherd et al. (1997), lie inside this clump. Nevertheless, we recall that this is a 
lower limit. Therefore the possibility of rotation is as plausible as that of source multiplicity, and 
the discrimination between them should perhaps be left to further high-spatial resolution images of 
either C'^O or high density tracers such as CH3CN. 
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Fig. 2. Molecular outflows: the solid lines represent the blue wings and the dashed lines the red 
wings of the '^CO(2-l) emission, with contours starting from 5cr and increasing by steps of 2<t 
(G43), 3cr (G10.47, G16, G19, G28, G48, G75) and 5ct (G10.62, G35). The rms values of flie blue 
emission are 0.2 (G10.47, G35), 0.3 (G10.62, G16, G19) and 04 (G43, G48, G75) K. The rms 
values of the red emission are 0.2 (G10.47, G16, G35), 0.3 (G10.62, G19, G28), 0.4 (G48) and 
0.5 (G43, G75) K. The blue-wing and red-wing velocity ranges used for integration are shown in 
TablelD The grey scale depicts the correspondingC'^0(2-l) emission from 5crby steps of 5cr. The 
rms values of the €'^0(2-1) emission ai-e 0.2 (G10.47, G19), 0.3 (G16, G35), 0.4 (G10.62, G28, 
G43) and 0.5 (G48, G75) K. The C"^0(2-l) half -power contour is drawn in white. Triangles and 
squares represent CH3OH and H2O maser spots, respectively (from Walsh et al. 119981 Hofner & 
Church well. 1 1 9961 and Forster & Caswell, ! 19991 1. and crosses mark the position of UC Hii regions 
(from the authors indicated in column 6 of Table [T]|. The axes show offsets in arcsec from the 
absolute coordinates given in Table[T] We note that the blue emission peak indicated with an arrow 
as B in the man of G48 corresoonds to a difi'erent clumo with a slishtlv bluer svstemic velocitv 
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Fig. 3. Velocity gradients: the grey scale represents the velocity of the C O gas by steps of 
0.2 km s The bluest (white) and reddest (dark) levels correspond to 66.6 and 67.8 km s ' 
(G10.47), -3.35 to -2. 15 km s"' (G10.62), 33.45 and 34.75 km s"' (G35.20), 54.1 and 55.3 km s"' 
(G43.89), -0.7 and 0.7 (G75.78) km s '. The blue (solid) and red (dashed) contours of the 
'^CO(2-l) outflow maps are overlaid for comparison. The symbols have the same meaning as 
in Fig.m 
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Fig. 4. Top: Mechanical luminosity of the outflow versus bolometric luminosity. Middle: 
Mechanical force of the outflow versus bolometric luminosity. Bottom: Mass loss rate versus bolo- 
metric luminosity. Filled symbols correspond to the sources in our sample and open circles to data 
from Beuther et al. (I2002bl l. The vertical bars indicate, respectively, the variation of Lmec, ^^out and 
M with inclination: the lower limit corresponds to an inclination of the outflow axis of 30°with 
respect to the line of sight, whereas the upper limit corresponds to an inclination of 60°. For our 
data, the optical depth is taken into account (circles; see Tables |5] and |6]l. G28.87 and G43.89, for 
which no t,^ estimate has been possible, are marked as filled triangles and represent lower limits to 
the outflow parameter plotted. The dashed lines represent the best fits to the data compiled by Wu 
et al. ( |2004||2005"l l. Note that these are not fits to the data plotted in this figure. 
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Fig. 5. Top: Kinetic energy of the outflow versus bolometric luminosity. Middle: Total outflow 
momentum of the outflow versus bolometric luminosity. Bottom: Outflow mass versus bolometric 
luminosity. Filled symbols correspond to the sources in our sample and open circles to data from 
Beuther et al. (I2002bl l. The vertical bars indicate, respectively, the variation of Ekin and p with 
inclination: it corresponds to an inclination of 60°with respect to the line of sight. For our data, 
the optical depth is taken into account (circles; see Tables |5] and |6]l. G28.87 and G43.89, for which 
no Tn, estimate has been possible, are marked as filled triangles and represent lower limits to the 
outflow parameter plotted. 
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Fig. 6. Ionising photon rate, A^Ly against bolometric luminosity, Lboi- Filled circles correspond to 
our data and open circles to data from Beuther et al. (I2002bl l and Sridharan et al. (I2002I I. The dashed 
line represents the expected curve for an individual source on the ZAMS (Panagia [1973l l. The arrow 
indicates the direction along which the points shift if the distance to the source is modified. 




Fig. 7. Mechanical force of the outflow versus Lyman photon rate of the UCHll region(s) associ- 
ated. Filled circles (corrected for Tm) and triangles (with no Tm estimate) correspond to the sources 
in our sample and open circles to data from Beuther et al. ( |2002b| ) and Sridharan et al. ( 120021 ). 
The vertical bars indicate the variation of with inclination: the lower limit corresponds to an 
inclination of the outflow axis of 30° with respect to the line of sight, whereas the upper limit 
corresponds to an inclination of 60°. 
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Fig. 8. Top panels: Total molecular outflow mass versus clump mass from our C'^0(2- 1) data (left) 
and from continuum (sub)millimetric surveys by the authors cited in Table |2] (right). The dashed 
line is the best fit obtained by Beuther et al ( I2002bl ) for their data. Bottom panels: Mechanical 
force of the outflow versus core mass as derived from the C'^0(2-l) emission (left) and from 
continuum (sub)millimetric surveys by other authors (right). The vertical bars indicate the variation 
of Font with inclination: the lower limit corresponds to an inclination of the outflow axis of 30°with 
respect to the line of sight, whereas the upper limit corresponds to an inclination of 60°. Filled 
circles represent the sources in our sample and open circles data from Beuther et al. (2002a), for 
which clump masses correspond to data from continuum emission at 1 .2 mm in all four panels. For 
our data, the optical depth is taken into account (circles; see Tables |5] and |6]l. G28.87 and G43.89, 
for which no estimate has been possible, are marked as filled triangles and represent lower limits 
to the outflow parameter values. 
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Fig. 9. Mout/Mciump versus iboi/A/ciump- Filled circles correspond to the sources in our sample (only 
those with a Tn, estimate are represented) and open circles to data from Beuther et al. (I2002bl l. 
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Fig.A.l. Comparison between the '^CO(2-l) and C"*0(2-l) spectra from the central core in 
the map of G 10.47 (upper panel), and the corresponding spectra seen towards the secondary red- 
shifted core in the same map (lower panel; G10.46B core), located at an approximate offset of 
(-30", -25")- The vertical dashed lines indicate the respective Vlsr listed in Table [3] Note the 
difference in their respective systemic velocities for both the '^CO(2-l) and C'**0(2-l) lines. 



